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ABSTRACT 


Cuts were made in a 840 mm x 25.4 mm x 5 mm Polymethyle- 
methacry late (P.M.M.A.) beam which were joined using Ethylenedi- 
chloride. The joint thus formed is ultrathin and can be 
approximated to be of zero thickness. Viscoelastic property of 
interfacial region has been varied by adding plasticizer Dioctyle 
Phthalate to Ethylenedichloride. For controlled force applied, 
transverse wave's acceleration (measured at different locations on 
the beam) consists of the forward propagating wave, reflected wave 
from joint, reflected wave from end and doubly transmitted through 
the joint; and the reflected wave from other end. At any given 
location these waves are seperated in time and Fast Fourier 
transform of these waves was taken for locations before and after 
the joint. Results show that the amplitude gain in dominant 
frequency decreases as the plasticizer content in the joint is 
increased. By mechanical testing of specimens made of similar 
joints it has been found that ultimate strength of the joint also 
decreases with plasticizer content. 
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CHAPTER 1 


Introduct ion 

The use of adhesive bonds in primary structures - is restricted 
due to limited growth of methodologies for predicting their life. 
It is desired to develop nondestructive testing (NDT) 

to establish quantitative relations between the materials 
structure and its mechanical performance . Among various NDTs. 
the commonly used techniques are pulse echo through transmission, 
C scan, resonance, spectroscopy, vibrations under impact etc. All 
nondestructive testing procedures proceed in two stages : in the 
first one (sometimes called the forward problem) we determine wave 
speed for given joint properties. In the reverse problem, we make 
measurements of wave speed through joints and use this information 
to find out the joint properties j( i . e . , the adherend thickness and 
wave speed in it, the adhesive thickness and the wave speed in it 
and the modulus). 

The important characteristics of joints is that there is a 

definite limit to the bond strength developed between specified 

6 1 0 

adherends for a particular adhesive . If it is decided to use 
double overlap joints, the load capacity increases with increase 
in overlap until a \ior>j'Vin^ value is attained. Beyond this overlap, 

no- .. load transfer can be achieved. Experiments have shown 

* 

that load transfer is mostly confined to the two end zones with a 
1 ight ly- 1 oaded elastic troughs in between. The second inherent 
characteristics of the joint is that the maximum potential bond 
shear strength is determined by the adhesive strain energy in 
shear per unit bond area. This is shown to be very important 
criteria. The precise shape of the stress-strain curve has been 
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shown to have little effect on the lie* tine joint strength. Joint 
efficiency charts are provided for common adherend materials. 
Since the adhesive films behave differently under peal and shear, 
the results shown in them can at times be misleading. For moderate 
thicknesses, the limiting shear strength of the adhesive may limit 
the joint strength. As opposed to this, for thicker adherends, the 
limiting strength is invariably the peal strength. 

When a joint undergoes a fracture, the failure could occur 

either because the adhesive has failed or the interfaces between 

11-21 

the adhesive and the adhearend has failed . The former is 

called cohesion failure to distinguish from the latter called 
adhesion failure. In an actual joint the failure normally occurs 
by a combination of both these. The failed joints have been 
examined under the microscope and invariably the failure is 
initiated at the microvoids created either by the air entrapped or 
the reaction products released during the curing of the adhesive^ ^ 

Commonly encountered bonding problems can be classified into 

following types : debonding, cohesive weakness and interface 

weakness. Debonding includes voids and complete sepration of the 

adhesive and adherends. This situation can normally be detected by 

the techniques of pulse echo, through transmission, C-scans etc., 

cohesive weakness by techniques aj resonance, spectroscopy etc. As 

opposed to these, the detection of interface weakness is more 

difficult. The adhesive and adherends could appear in physically 

perfect bonding state, but could lack the strength. Traditional 

ultrasonic techniques do not distinguish between good and bad 

bonds. Work has been initiated in the development of new 

techniques and theoretical models to understand the wave 

. 28-30 

propagation at the bond interface. Various theoretical models can 
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be classified into two types 25 29 t (a) a complete bond thickness 
weakness model and (b) interface wes‘;r:iss model. 

The bond thickness weakness model refers to wave propagation 
(e.g., plate waves, leaky waves, interface waves etc.) through the 
adhesive layer and two larger adherends. The received signals are 
affected both by interface as well as by cohesive weakness. For 
the second type models, the interface weakness is evaluated more 
directly by examination of the transmission factors from the 
interface between the adhesive and the adherends. The existence of 
an interface wave has been demonstrated which has smooth contact 
with the adherends and the adhesive and acts as a viscous liquid 
layer. In another quasi-static model of a weak bond, a spring with 
an elastic constant and mass was assumed to connect two solid semi 
infinite spaces. In these models, different bonding qualities 
refer to different bonding conditions, such as (1) debonding, (2) 
kissing bonding (3) smooth (or slip) bonding (4) general weakness 
bonding and (5) rigid (welded or perfect) bonding. 

The rigid bonding is considered as a perfect contact between 
two media and it assumes the continuity of traction force and 
particle velocity at the interface. The kissing bond assumes that 
the traction force and particle velocity are continuous only in 
the direction normal to the interface and that the tangential 
components of the traction force is zero. The reflection and 
refraction factors from a kissing bonds come out to be of 
frequency dependent mathematically which unfortunately is contrary 
to the experimental finding. The smooth bonding is of primary 
concern because it can easily be missed with a traditional normal 
beam longitudinal wave inspection procedure. The smooth bond can 
be easily simulated experimentally with a thin film between two 
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solid media. A perfect bonding in an adhesive joint can never 
exist and some of the major parameters that affect the bond 
strength are bondline thickness, surface roughness, chemical 
conditions, adhesive type, environmental conditions and the aging 
effects. The boundary layers have different physical and chemical 
properties . 

The fracture mechanics approach has been widely used to 
correlate the quasi-static and fatigue crack growth behaviour of 
the adhesive joints. The fracture parameters are often calculated 
using finite element methods because the governing equations for 
arbitrary joint configuration loading modes in complex structure 
cannot be solved analytically . One way is to treat all planar 
adhesive joints as comprising of one or more sandwitch elements 
and the free body diagram of each written. The reactions in each 
adherends are determined and stresses in each adhesive layers 
obtained using plate bending analysis. 

In this work, we propose to use flexural waves generated 
through impacting for characterizing the joint interface. In doing 
this we have made the following two kinds of measurements of the 
flexural waves. 

(a) The measurement of acceleration at a given point. In this 
process we have argued that the [forward and reflected transverse 
waves are separated in time. 

(b) Comparision of forward waves at immediate points across 
the joint. We have taken a poly methyl methacrylate (PMMA) beams 
of uniform cross section and produced a cut on it. The two pieces 
were joined using Ethylene dichloride solvent. The viscoelastic 
nature of the joint was varied using & suitable plasticizer. 
Through our experimental analysis we have developed a correlation 
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between the joint property and transverse waves generated in it. 


Further, we have also solved problem of a viscoelastic beam 
clamped in middle and vibrated by a sinusoidal force^^ , solution 
has been used for material characterization. 




- - CHAPTER 2 


SOLUTION 0£ FLEXURAL WAVE PROPAGATION IN VISCOELASTIC BEAMS A 

FOURIER TRANSFORM APPROACH ♦ 

Let us consider a uniform viscoelastic beam with a simusoidal 
force applied at its center. Because of the symmetry, we solve the 
problem of beam with force applied at one end as shown in figure 
Ar. The governing differential equation for transverse displacement 

y is pyurutnr+l ^ 


dV _ dhk a 2 y 


6X 


ax' 


= F ( x ) - m 


at 


(i) 


For the viscoelastic material, the stress and strains are related 
in general by (_ *1 


[ Q o + 37 + «* < ^ > * '■■■■ S. J,] *(t) 

= [ p o + p ! a?*- S] E<t) 


( 2 ) 


where a(t) and e(t) are the stress strains at any given and 

ovot^rnoi 

to Sn and P o to P n known ^constants . We define o and e in 
Fourier transform domains as o (o) and £ (o ) as follows 


oo 

a(o) = J o(t) e -i 

—CD 


-iwt 


dt 


(3) 


£<«) = f e(t) e" iWt 


dt 


(4) 


we 


obtain j~rcxm 
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5(o) 


E (o).e(t) 


(a) 


where 


r *, u \ - p < i0) > 
E <U) ‘ Q(t») 


(b) 


P( io) 


p Q + Pj(io) + r^rrcr)-' 


. P ( io) 
n 


n 


(c) 


Q( io) = q q + Q 1 (io) + Q 2 (io) 2 + Q n (io) n (5) 

One can show that in this domain, eqn(l) reduces to 


E* I f m (io) 2 y = F(o) (6) 

dx 4 

where F(o) represents the Fourier transform of the -|-o tcjz . 
appl i ed . 

Now we consider the problem of a beam held in a clamp of mass 

M. The clamp is supported in such a manner that it is able to move 

freely in the y direction but can not rotate in the x-y plane. It 
is assumed that a sinusoidal force acting upon the 

clamp, as shown in figure 1. The boundary conditions for this 
systems are(a) zero bending moment at x = 1 (b) zero shearing 

force at x = 1 (c) zero slope of displacement at x = 0. The last 

boundary condition is the equality of shear force at x = 0 to the 

inertial resistance of clamp. 

In the transform domain these can be written as at x = 1 


Moment M = E I 


2 - dx2 

*4 = o 

dx 


= 0 gives 


(7) 
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Shear force V = 


* d^y 

E I ¥r = 0 gives 

dx” 5 


At x = 0 


Slope = 


,3- 

^ = 0 
dx J 


¥ = o 

dx 


( 8 ) 


(9) 


Similarly, the equality of shear force and the inertial resistance 
of clamp gives 


■M$)- 


-■($) 


= M( iw) ^y 


dt 


( 10 ) 


■where M is the mass of clamp. The differential equn. in the 
transform domain becomes 


dx 


*-[£]* ■ i - 


,® F e 6(x-0)e * <l) * dt 


E* I 


or 


.4- 
— <l y 

dx 4 


a 

<N 

3 


r zic i 

L „* 

y = 

1 * 

L E I J 


L E I J 


S(x-O) 6(u-« ) 


( 11 ) 


We define k as 


.4 _ m<o 
: “ * 

E I 


( 12 ) 




[ 1 5(x-0) «(«-« ) (13) 

L E I J ° 

The boundary condition in eqn. 10 is rewritten in terms of k as 
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<*LZ = 

dx 3 


— I k^y at x = 0 

m J 


( 14 ) 


The general homogeneous eqn. corresponding to the problem of eqn. 


13 is 


e*i d *y - k 4 ; 

dx 4 n " 


whose solution is 


y =A Sink x + B Cosk >< + F Si nhk + G Coshk ( b ) 

n nnnnna nn 


Now, we apply the Boundary Conditions 7 to 10 to obtain 


(Cosa + Cosha) A + (Sina) B + Sinha G = 0 

n n n 


(17) 


(Sina + Sinha) A^ + Cosa - (Cosha) G^ = 0 


(18) 


««A n + G n 


(19) 


where P - 

ml 


a = k 1 


( 20 ) 


For getting non trivial solution for A , B and G we must have 

n n n 


- (Cosa + Cosha) (Sina) Sinha | 


(Sina + Sinha) . Cosa - (Cosha) 


( 21 ) 
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This sets a condition on a and 0 such that 


» «0 


( S ina ) (Cosha ) + Cosa Sinha 
1 + Cosa Cosha 


( 22 ) 


This equn. has theoretically infinite solution for a for given 0. 
Each value of gives a solution y n » which after inversion will 
give the nth mode shape y ft (x) in eqpi. 15. This can be derived as a 
viscoelastic beam of above type. 

y n (x) = A n ( Sink^x - Sinhk n x) + B^ Cosk^x + G^Cosh^x 

(23) 


where 


(Cosa + Cosha ) 
n n 


n 


B 


^(Sina^ + Sinha^) - ( 2/a^0 )Cosa^J n 

(Cosha + (a 0/2) (Sina +Sinha )1 
L n n n n J 


n 


I a ft _______ 

Cosa - . n x (Sina + Sinha ) 

L n (— =— ) n n J 


n 


(a) 


(b) 

(24) 


Above G_ i§ unspecified, but is used to normalize y (x). Hence we 

n 

note that one solution of equn. 22 is a = 0. But if a is equal 

o o 

to zero, this would correspond to k = 0 and 
y o (x) = B q + G o 

which is constant. This represents a linear translational of 
motion of the beam in phase with the clamp (no bending). 


In order to determine Q n in'eqn. 15, we notice that in the 
transform domain, it can be written as 


oo __ 

y (x,w) = £ y„(x) Q t> < ( °) 

n=0 n n 


(25) 


If this is the complete solution, then it must satisfy the 
PDE in eqn. 13. Or 


m .4- 
® d y _ 


CD 


2«F 6(o)-w ) S(x-O) 


n=0 dx 


T % - X 4 Z - ° ■ ° 


n=0 


n n 


* 

E I 


(26) 


Combinin this with the homogeneous equation in eqn. 16a, we get 


CD 


. . 2 «F Sfo-u ) 6(x-0) 

E <*i - * 4 > ?» % = ° . 

n=0 El 


(27) 


Multiplying this by y^(x) and integrating from x = 0 to x = 1 we 

get 


00 , A _ r l 2?tF y.(o)6(«-o ) 

E <k 4 - k 4 > yy .dx = -2-^ 

n — » J Q J r 

n-U o El 


(28) 


We know that the shapes of the modes y and y. are such that 

n j 


J y n (x) yj(x) dx = - (M/m) y n (°)yj(°) 

n * j 

Substituting these in eqn. , we get _ 

- - X 4 ) *n<°> * <- 4 - ^ 4 ) 


(29) 


T 


n t j 


2nF m 
o 


( 30 ) 


_ — 6 (o - o ) 

ME I ° 

The relation in eqn.(30) is actually an infinite set of linear 
simultaneous equations as follows 


. Q mk . . _ 

" k * M yo ~ * k l “ k ) Qi y l<°) 


- k4) Q n y n<°> 


2nF m 
o 

ME* I 


6(0 - o ) 


(a) 


A a i mk 

(o - k 4 ) + (k 4 - k 4 ) — - 

Vo , VK 1 K ' M y. (o) 


(k 4 - k 4 ) Q y (o) 


2nF m 
o 

ME* I 


6 ( 0 ) - 0 ) ) 
o 


(b) 


(O - k 4 ) Q o y o - (kj - k 4 ) Q 4 y A (o) 


. . Q mk 

(k 4 - k 4 ) j - T 

n My (o) 


n 


2nF m 

— j 5 - 6 (o - o ) 

ME I ° 


(c) 


Substituting the nth eqn. from eqn. ( 1 ) gives the following 


(o - k 4 ) y Q [ (mk/My^) + 1 } 

(k i ~ k4) y n (o) 1 ® k /My„(o)] + 1 


(31 )_ 


Now we determine Q q since zeroth mode represents the rigidbody 
motion, this simply means 
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F e 
o 


i<o t 
o 


= (M + ml) d y( V J - 


dt 


(32) 


In transform domain this becomes 


2nF Q 6((J - Q o ) = (M + ml ) ( iw) 2 _ = M [l + gp] 


(o - «‘> y Q Q c 


(33) 


giving Q q as 


Q, 


2nF S(o> - « ) 
o o 


M 


(ml /M) +1 y Q (o - w ) 


(34) 


On combining eqns . (31) and (34) and inverting gives 


Q (t) 
n 


* F e‘“»‘ 
o 


M(y (o) |(o p +iw ni ) 2 - 

n L ni o 


U ♦ -®| — ] 

L My Z (o) J 
n 


-1 


This means that the transverse deflection is given by 

iw t 




y(x, t) = £ 


y (xTf e ° 


n o 


n=0 M(y 


(o) 

n 


(to + i<o . ) ' 
nr ni 


o 


[i ♦ -a|— ] 1 

L My‘(o) J 


(35) 


f y^x^-dx = k 
j n, 


We find the accelaration from deflection by differentiating 
it y twice with respect to time as 
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o 


CD 

a( x , t ) = - l 

n = 0 


y ( x ) F w e 
n o o 


r 2 .8 

I « 


2 ..2 


,8 - 

1 ta 


M y n (o) li j (Ej + iE 2 ) -« G 

[ 


x 1 + 


A(x) 


My‘(o) 


( 37 ) 


y (x) F w 
n o o 


L 1 ai J L o My ( o ) J 

n 




B( x) 


f i!a! 

t - « 8 1 

n 


iii ) 


2 ^2 

— (j - e + F 
2 ”o ^1 b 2 


) * «i 


2 e 2 

2 




od f 8 2 ^ 

= J o B(x) [ u 2 - E 2 ♦ E 2 j Cos^t-lE^Sin^t] 

n 


r ® _ r i 8 2 

+ i E B(x) -Ls- 

L n-0 L a* I 


n 


+ 2E i E 0 CosO) t 

I L O 


]] 


2 “o - E l + E 2 ) S,1 >VJ 


(39) 


If JLf then the imagenary part of eqn . (39) will give 

the a(x,t) as 


CO 


a ( X , t ) 


= E "B(x) [ -■■■?, , u 2 - + eI 1 Sinu t| 

n=0 L a 8 \ L ° 1 2 J ° J 


n 


+ 2E 1 E_ CosM t 
XL 0 


]] 


(40) 
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( 41 ) 


a(x, t ) = E B (x) Sin(ot+$) 

n=0 n 


i(x,t) = Amp Sin(ut+<f>) 


where tan<#> 


2E 1 E 2 


( 77 Jil - A * 4 ) 

n 


CD 

Amplitude = E ® (x) 

n 


( 42 ) 

( 43 ) 


( 44 ) 


Now if we know the Amplitude of the accelaration at any two 
distance and from the clamsj then we can get values of Ej & 
E, corresponding to the frequency « . 
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CHAPTER 3 


EXPERIMENTAL SEIUP AND PROCEDURE FCR JOINT 

CHARACTERIZATION 


A P.M.M.A. rectangular beam (length 840 mm, width 25.4 mm 
thickness 5 mm, made with two pieces (300 x 25.4 mm x 5 mm and 540 
mm x 25.4 mm x 5 mm) and joined together with diethylene chloride 
is hung from a fixed iron frame with two thread* tied to the ends 
of the beam. The points of contacts between the thread and the 
beam and between the thread and the iron frame are rendered 
immobile with an epoxy glue so as to obviate any friction which 
might interfere with the ideal propagation of flexural waves in 
the beam. To generate flexural waves an impact hammer is hung from 
the iron frame with two threads forming a V shape, the motion of 
which is re»4r~icted to lateral dimensions only, with respect to 
the beam. Quartz piezotron accelerometers serve to pickup the 
acceleration response at different points on the beam located on 
either side of the joint. A multichannel signal coupler conditions 
the accelerometer response to a degree commensurate with the two 
F.F.T. nalyzers. A movable support serves to enhance the 
reproducibility of the impulse force generated by the impact 
hammer when released from a given height. A small cube of 
(P.M.M.A.) negligible mass affixed in the middle of beam at the 
point of impact to make impulse force spatially localized to that 
point only. An I.B.M. compatible personal Computer (PC-XT) accepts 
the sampled data transferred from the F.F.T. analyzer through the 
General Purpose Interface Bus (G.P.I.B.), the pertinent details of 
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BLOCK DIAGRAM OF EXPERIMENTAL SET UP 






FIG 2 EXPERIMENTAL SET UP FDR JOINT CHARACTERAZATION 
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FIG 3 SIDE VIEW OF EXPERIMENTAL SET UP 


FOR JOINT CHARACTERIZATION 



which are provided in Appendix* An application program in Basic 
interprets the sampled F*F*T. data which is then processed with 
the P.C*D.S*P« (Personal Computer Digital Signal Processing) 
software* app 1 i ca t i on program 

is given in Appendix^ • 

Accelerometer responses at two symmetrical points with 

respect to impact position are choosen such that one of them lies 
just before the joint, other one lies symmetrically towards the 

opposite end* When these two positions responses were compared, 
we found that after about 0*9 miliseconds their nature changes* 
This shows that after 0*9 miliseconds of impact the point which is 
just before joint starts receiving reflected wave from joint, or 
in 0*9 miliseconds forward wave has travelled from point of impact 
to the point before joint* Since for about 1 miliseconds the wave 
is going only in forward direction, therefore we have left data 
for first 1 milisecond and we have taken F. F*T. of next 4 

milisecond data only* W c have found that this zone contains 

* • • k 
precious information about the joint char act eres"tt<^|^ call this 

n. A 

zone as pristi* e zone. 

A 

Procedure 


(a) The acceleration responses as well as the impulse force 
response of the impact hammer at points mentioned in Fig.l are 
taken for different joint properties and different hammers. 


(b) With the hammer tips of varying softness and for different 
joints the acceleration response^ at two points one being 5 cm 
before joint and another 5 eras, after the joint are stored in 
F.F.T. analyzer and analyzed. It has been found that lower 
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frequencies in succession get amplified after the joint while the 
higher ones in succession get accentuated. For beams with joints 
of different D.O.P. percentages the anpd i tudo- (pa-in that is ratio 
of amplitude of same frequency before and after the joint were 
noted and it was found that amplitude gain across the joint of 
major frequency is a strong monotoni cal ly decreasing function of 
D.O.P. percentage. The major frequency is found to be same before 
and after the joint. For higher frequencies amplitude gain is more 
or less one this shows that at higher frequencies joint being very 
this acts 1 Hre a solid beam material itself. 
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CHAPTER 4 

EXPERIMENTAL SETUP AND PROCEDURE TO CHARACTERIZE 

THE MATERIAL 


The polymer beam shown in Fig. 2 is clamped in the middle and 
clamping device is connected to the strain gauge based force 
transducer mounted on the plate, which is connected to shaker 
head. Force transducer senses force in terms of strains, b^frng 
caused in strain rossets, mounted in force transducer. The desired 
frequency is generated by a function generator whose amplitude is 
amplified by a power amplifier which drives an electromagnetic 
shaker. The force to which the beam and clamping assembly is 
subjected is sensed by the load cell (force transducer). The 
strains caused due to force change the resistance of the strain 
rossets fitted inside the force transducer. Dynamic force was 
measured using the force transducer connected to a Q-unit of 
oscilloscope. 

Two piezotron quartz accelerometers are stuck to the beam at 
distances 15 cm and 30 cm from the clamp, which senses the 
accelerations at the corresponding points and send the signal to 
F.F.T. analyzer which displays the time data on the screen. The 
corresponding acceleration and force from the given voltage on the 
screen is found by multiplying it with the calibration factors 
given the corresponding manuals. 

Procedure 

Desired frequency at which force has to be applied is decided 
and then that frequency was generated by a function generator and 
signal was amplified by power amplifier. Then that signal was send 
to electromagnetic shaker to vibrate it with that frequency. Force 
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EXPERIMENTAL SET UP FOR MATERIAL CHARACTERIZATION 






applied by shaker on load cell as well as clamping device, with 

beam clamped in it is sensed by load cell and displayed on 

osciloscope screen as a function of time in terms of milivolts. 

Since according to analysis shown in chapter 2 fpr finding out 

modulus of elasticity^ only amplitude of force is desired 

therefore we read amplitude of force from screen of osciloscope in 

terms of milivolts and then this value is multiplied by 

corresponding calibration factor that is 12*83 N euitons to get 

mv 

force in M-au^crn^ * Similarly! accelerometer response was also 


di spl ayed 

on 

F.F.T. analyzers screen in 

terms 

of milivolts. 

Ampl i tude 

of 

acceleration response were 

read 

f rom screen i n 

milivolts 

and 

converted in acceleration 

units 

2 

(m/ sec ) . Now 


knowing amplitude of force and accelerations at two points and 
mass of clamp and mass of beam, length, width and size of beam 
etc. and using the expressions developed in chapter 2 we can find 
out the complex modulus for given frequency by solving those two 
nonlinear expressions by iteration. Similar procedure was repeated 
for different frequencies and corresponding complex modulii were 
found. Since the equations were simultaneous nonlinear equations 
they resulted in many solutions. Out of those one solution which 
was reasonable was 

Ej = 5.87 x 10 10 N/m 2 
E 2 = 1.12 x 10 8 N/m 2 

This value remained more or less same for higher frequencies say 
100, 200 etc. Further higher frequencies could not be generated 
with our instrument we were not able to get E* vs frequency 
relation for P.M.M.A. beam only in frequency range 50 Hz to 250 Hz 
we were able to find its values. 
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CHAPTER 5 


EXPERIMENTAL SETUP Al© PROCEDURE FOR ESTIMATING 
ULTIMATE JOINT STRENGTH 

An Universal Testing Machine (U.T.M.) stretches the P.M.M.A. 
specimen as shown in Fig. 3. Specimen is made as per A.S.T.M. 
standards as shown in Fig. 4. Specimens were stretched at slowest 
possible rate that is at 0.2 mm/min. because specimen is made of 
viscoelastic material whose modulous of elasticity depends upon 
rate of loading. Its ultimate strength also depends upon rate of 
loading. Minimum speed was selected so as to get highest strength. 
Specimens were stretched till failure. Simultaneously it plots 
load vs extension curve. Salient features of Testing were 
Rate of stretch = 0.2 mm/min. 

or cross head speed. 

Chart speed = 1 mm/mirj. 
load capacity of machine = 5000 kg 
specimens were 15 mm wide and 5 mm thick. Because jaw of machine 
cannot accomodate width more than 25 mm and thickness more than 
12.5 mm. 

Procedure 


For specimens with different D.O.P. percentage stress taken 
at the time of failure was computed and averaged over fine samples 
of same specimen. The ultimate strength is found to be a strong 
monotoni cal ly decreasing function of D.O.P. percentage. 


2 * 


Experimental Data 

at <J =60 Hz (P.A. = Peak Amplitude) 


at 

. 15 



P.A. 

2 

327.544 mV 

at 

.30 



P.A. 

= 

709.679 mV 



at 

0 ) 

0 

= 30 Hz 



at 

.15 



P.A. 


3.25906 V 

at 

.30 



P.A. 

= 

1.04814 V 



at 

(0 

0 

= 300 Hz 



at 

.15 



P.A. 

= 

1.33201 V 

at 

.30 



P.A. 

= 

600.498 mV 



at 

(0 

0 

= 100 Hz 



at 

. 15 



P.A. 

= 

2.38561 V 

at 

.30 



P.A. 

= 

1.5886 V 


f ■ 

at 

(*) 

- 0 

= 200 Hz 



at 

.15 



P.A. 

s 

595.039 mV 

at 

.30 



P.A. 

= 

185.609 mV 


Computational work for calculation of E^ & E 2 
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CHAPTER 6 


Results and Discussion 

The Joints in a primary structure could be mechanical welded 
or adhesive in nature. In this work we focus our attention on the 
adhesive joints in polymer beams. The polymer consist of large 
molecules which may have linear or cross linked chains. The latter 
are network in structure where all polymer molecules are linked 
together giving a giant molecule. The difference between the 
linear and network molecules lies in their solubility in that the 
latter does not dissolve but incorporate solvent into themselves 
through swelling. 

When we wish to prepare joints in polymeric beams, we have 
only the following limited options. 

1. Use a suitable adhesive which sticks the two adherends. 

2. Use a solvent system which disolves the material of 
adherends and its evaporation leaves behind a joint. 

In any given polymer structure (as shown in figure 4) the chain 
ends of polymer molecule are dangling out when we apply an 
adhesive the polymer chain ends must unentangle and penetrate into 
the adhesive and the driving force for this is obtained by the 
favourable difference in their chemical potentials. This 
inaterpenetrat ion gives rise to two interfaces as shown. In 
contrast to this, when a solvent (is used the chain ends of the two 
adherends must dissolve into the solvent and interpenetrate into 
each other. On evaporation of the solvent, this leads to single 
interface. The strength of the joint arises because of this 
interpenetration of chain endseither into each other or into the 
adhesive. The complex combination of the strength of the adhesive 




and that of the interpenetration of the polymer chains is not yet 
understood very well and the purpose of this study is to throw 
some light upon this. 

In this study we have chosen a beam made of 
Polymethylmethacrylate (PMMA) and we have cut it at a suitable 
location. The two pieces of the adherends are joined using 
Ethyl enedi chi or ide solvent. This has very high vapour pressure and 
gives an excellent joint. In order to vary the property of this 
joint we have mixed Dioctylphthlate (DOP), a plasticizer, (CgH^ - 
1,2 - { C0 2 CH 2 CH (C 2 H 5 ) (CH 2 ) 3 CH 3 ) 2 with it while making the 
joint. The result of this way of joining is that the solvent 
evaporates with DOP remaining within a very thin layer around the 
interface. It is desired now to determine the changes of the 
viscoelastic properties of this interface. 

For determining the properties of the interface we have 
employed the use of transverse flexural waves generated through 
impacting. We strike a beam with a small metallic bob by swinging 
it from a known height it generates flexural waves on impacting 
which travels both ways with equal likelihood as shown in Fig. 5. 
Let us say we decide to measure acceleration at point A using a 
suitable accelerometer. The forward wave from the point of impact 
travels to this point first. After it has travelled across and 
reached point of joint J, some of it is reflected and the 
remaineder is transmitted through it. This is reflected from the 
end of the beam B passes through the joint once again and 
ultimately reaches the point of measurement. Similarly energy 
travelling towards C is also reflected in due course and reaches 
point A. Since the waves take a finite amount of time to travel, 
it is expected that these reach point A at different time and this 
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way produces a separation between these. This separation is not 
expected to be clearcut, since higher frequency waves travel 
faster and may overlap with other waves. However, it appears 
possible to distinguish the various regions. 

Figure 7 shows the pristive time zone. The F.F.T. of which 
carries the desire information about the viscoelastic properties 
of the interface which acts as a natural band limited filter 
modifying the elemental frequencies. 

Figure 7 shows acceleration response before joint in time 
domain and Figure 8 shows acceleration response after joint in 
time domain. Figure 9 and Figure 10 show F.F.T. of acceleration 
response before and after joint. Figure 11 shows amplitude of 
dominant frequency vs different locations on the beam. 

A drop of eighty percentage in amplitude gain of the major 
frequency of 2517 radians per second, shown in figure 12, when 
D.O.P. percentage in the joint is varied from 0 to 20, reflects 
strong relation between the amplitude gain and the viscoelasticity 
of the joint. 

With linear regression the following first order model fits 
approximately, the experimental data. 

y 1 = 25 - x (45) 

where 

y^ = amplitude gain; 

x = % DOP. 

The relationship between ultimate strength and percentage 
D.O.P. is modelled with the following first order polynomial 
approximately by 

y 2 = 250 - 10 x (46) 
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Fig 7 Acceleration Response 
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FIG 13 ULTJ 



yg = ultimate strength of joint; 
x = % DOP. 

Eliminating percentage D.O.P. we get, by solving the above 
two equations, the following relationship between ultimate 
strength of joint and the amplitude gain of the major frequency of 
the pristive zone. 

y g = 10 y 1 (approximately) (47) 

where 

y^ = ultimate strength; 
y^ = amplitude gain 


40 



mmm 


CHAPTER 7 
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■! 

Conclusions 1 


c*c{\u-c \vjz 

We observe that the strength of <<>.. joint consists of a 

A 

complex combination of the strength of the adhesive layer and that 
of the two interfacial regions. It is possible to study the 
strength of the adhesive material independently, but very few 
studies are reported concerning the interfacial region of the 
joint. We have chosen to study a beam of polymethyl methacrylate 
(PMMA) and joints made in it by using ethylene dichloride as a 
solvent. The thickness of this joint is negligible and can be 
assumed to consist of only interfacial region. We have chosen to 
vary its property by using dioctyl phthalate (DOP) plasticizer and 
wish to study its effect upon the flexural wave propagation 
through it. 

We observe the following : 

1. The flexural wave generated across the joint through 
impacting carries information about the joint while getting 
t ransmi 1 1 ed . 



2. Joint interface acts as a natural band limited filter. 

3. The frequencies across the joint are preserved while their 
amplitude gets modified. 

4. The amplitude of the major frequency of the prestive zone 
gets amplified. 

5. This amplitude gain is strongly correlated to the 


viscoelastic properly of the joint. 

6. The equation "^c-joyjcan be successfully used for the non 
destructive testing of the joint, fy - uJ-MinWie. V »y* 
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APPENDIX 1 

SPECTRAL ANALYSIS 


Several analytical techniques exist in the literature for the 
analysis of wave propagation problems. Central among these is the 
method of Fourier Transform analysis which is also sometimes 
called as spectral analysis. In this, any given signal is assumed 
to be formed by a superposition of many infinitely long wave 
trains of different frequencies. In the following we describe this 
in brief. 

CONTINUOUS FOURIER TRANSFORM 

The continuous fourier transform of a function F(t) is 
defined in the time domain -® to +oo as 



■fOO 

2 it F(t) = J C(u) e“ l0)t dw (b) (1.1) 

— CD 


Above €(<*>) is a continuous fourier transform, o> is the angular 
frequency and i is complex number equal to /~T. Eqn. 1(b) shows 

inverse transform and 1(a) shows forward transform. This arbitrary 

/ 

convention arises because the signal to be transformed are usually 
available in the time domain. The factor 2 n is for sequential 
recovery of forward as well as backward or inverse transform for 
the original signal. Using frequency relation 2itf = w, we get 


C(f ) 



-i2itf.t 

e 


dt 


(a) 
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F ( t ) 


df 


(b) (1.2) 


4-OD 

= J c* (f) e - l2nn 

~co 

where the superscript * refers to the fact that the variables are 
complex conjugate. Some of the properties of Fouriet Transforms 
are listed in Table Al . 

REAL TIME FUNCTIONS 

The following transform pairs are valid 


Fe( t ) ♦♦ 

1/2 

[ C(u) 

+ c(-tt) ) 

= Real [C(o>)] 

Fo( t ) * 

1/2 

l C(«) 

+ C(-o>) ] 

= i Image [C(w)J 


for both F(t-) and C(<<>) being complex. In practice the functions of 
usual interest in wave analysis are when F(t) is real. To see the 
effect of this, Eqns. are rewritten in terms of real (C R ) and 
imaginary parts (C^ ) of C in terms of F R and Fj of F as 


2n F r = J (C R Cos(ot) - Cj Sin (<*>t ) 1 dw (a) 

C R = J [F R Cos(o)t) - Fj Sin (ot) ) dt (b) 

2n Fj = J [C R Sin(wt) - Cj Con (ut ) ] dw (c) 

C R = J (F R Sin(ut) - Fj Cos (wt) J dN? (d) (1.3) 

Using the following decomposition 

Cos0 = (1/2) [ e l6 + e" 10 1 and (a) 

Sin6 = -i (1/2) ( e 10 + e -10 ] (b) (1.4) 
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If F( t ) is real only, Cj^ is even and Cj is odd mathemat ical ly 

C R ( “ U) = C R (W) ’ C I (" W > = "Cj (O) (1.5) 

which says that the functions are symmetrical and asymmetrical, 
respectively about the zero frequency point. This can also be 
expressed as saying that the negative frequency side of the 
transform is the complex conjugate of the positive side and is 
written as 

C (-u) = C* (-«) (1.6) 

FOURIER SERIES 

A fourier series representation of a function assumes the 
function as periodic in the time domain. Consequently, it can only 
be an approximation to the transient function found in wave 
analysis since these are assumed to be present on the full domain. 
In elementary terms, if the signal has large duration of zero 
amplitude and the analysis assumes it repeats itself on a period 
very large compared to the time of interest of the wave, then for 
all intents and purposes it behaves as an infinite wave. 


BASIC REPRESENTATION 

A function F(t) with period T can be expressed in the form of 
a Fourier series as 


F(t) 


OD 


1/2 a o + E 

n=l 


- 



r > 

m 

a Cos 
n 

/ 

H|r* 

v 

+ b Sin 
n 

2«n ^ 

t 4 



(1.7) 
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where the coefficient are obtained from 


a = | f F(t) Cos <2n n t/T) dt , n = 0,1,2, (a) 

n l 

T 

b = | f F(t) Sin (2n n t/T) dt, n = 1,2, (a) (1.8) 

n T o 

In these n can take on positive as well as negative values. It is 
apparent that a^ is symmetric while b^ is antisymmetric about n = 
0. Using the identities of Table A1 and introducing in terms of 
discrete frequency = 2nn/T, the above can be rewritten as, 


<0 


iu t 


F(t) = 1/2 a rt + 1/2 £ (» ~ ib ) e n 


o ’ tTi ' n n' 

n=l 


oo -i<D t 

+ 1/2 £ (a n + ib n ) e n 

n=l 


(1.9) 


T -iu t 

i + ib = £ F(t) e n dt, n = 0,1,2, (a) 

n n T J A 

o 


T t 

a - ib = £ f F(t) e n dt, n = 0,1,2,.... (b) (1.10) 

n n T J. 


Since the phase t is linear in n, the F(t) can be further 

■Written as 

-1 iu t ® i \ t 

F(t) = 1/2 £ (a n - ib ) e n + 1/2 a Q + £ U ~ i\) e 

n=® n " n=l 

( 1 . 11 ) 

This permits the simpler, more compact, representation of any 
given function, F(t) as, 

a 
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■Hx> i<*> t <t> icj t 

F(t) = 1/2 l U n - ib ) e n = l C e n 

—00 n 


—CD 


( 1 . 12 ) 


r T -iw t 

where C n = 1/2 ( a R ib n ) = 1/T J F(t) e n dt , 

o 


n = 0, ± 1 , ± 2 etc. 


from the above it is clear that 


431 *'■ 1 A " i 

-n n 


b = -b 
-n n 


n = 1,2,3, 


RELATIONSHIP TO THE CONTINUOUS TRANSFORMS 

Much of the literature on wave propagation use Fourier 
integral on continuous transforms. It is of interest, therefore to 
show the relationship between it and Fourier series. 

First, the series transform pair is rewritten 

CD iw t r T 

2* F(t) = l (T C ) e n Aw, T C = F(t) e n dt , 

_ n n j , 

-CD t 

(1.13) 

where w = 2rtn/T and Aw = 2*/T. Now let the time window encompass 
n 

from to consuming the signal is zero before time t = 0 ) . 

Thus Aw becomes dw and the summation becomes an integral, that is 

/ 

co iu t p +® — iw t 

F(t) + V T C cw e n dw, T C (w) * F(t) e dt 

n n J -® 

(1.14) 

These are the Fourier continuous transforms. In the region of 
applicability, they have the correspondence 
T C n 4 C(w) 

An alternative scheme to set up the correspondence is to 
consider the C Q component for pulse and how it relates to C(o) of 


cn 


the continuous transfora. 


DISCRETE FOURIER TRANSFORM 

The discrete coefficients in the fourier series are obtained 
by performing continuous integrations over the time period. The 
integrations are now replaced by summations as a further step in 
the numerical implementation of the continuous transform. 


APPROXIMATIONS FOR THE FOURIER COEFFICIENTS 

The essential starting point for the discrete transform is 
the series transform pair rearranged as 

® it) t T -i» t 

F(t) = | 1 C n e n , C n (w) = J F(t) e n dt 

—00 o 


Let the function F(t) be divided in M piecewise constant segments 

whose heights are f^ and base AT = T/M. The coefficients are now 
obtained from 


H-l +M/2 i». 

c „ ~ D „ = L f . J. - - « " = AT x '• 


= 0 


at/ 2 


m 


f Sin « n At/2 ) 


« AT/ 2 
n 


? e 


iw tm 
n 


= AT 


Sin « AT/2 
n 

w AT/*2 
n 


i f. 


- iw tm 
n 


(1.15) 


It is seen that this is the sum of the transforms of a series of 

at 

rectangles each shifted in times t = tm + j— • T * i% contribution 
of each of these will now be evaluated more easily. If n > m. 
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i.e., n = M + n f where n is an exponential terms because 


i<*> tn 
n 


tm 

n 

o 


- ioKO tm - in <t> tm 
e ° e ° 


( 1 . 16 ) 


Hence the s unmat i on contribution simply becomes 


M-l 

i f « e 
«=0 


-in w tm 


Looking at the other contribution, it is seen that the Sine 
function terms does depend on the value of n and is given by 


Sine (x) 


Sin(x ) 
x 


» X = « n &T/2 

= Kn f AT 
o 

_ _ 1 T _ „ n 

“ n TM' M 

Example 1 : 

Let us consider a superimposed wave having the following 
representation : 

y = aj Sin (w^t + ) + a 2 Sin (« 2 t + + a^ Sin (Wjt + $g)+b+z 


where, 

b 


a l ’ a 2 * *3 


V V *3 


constant representing d-c component 
amplitudes of waves , u 2 , « 3 frequencies in 
order <^< « 2 < 


s The phase differences. 

*■ * a random sriable representing noise (assumed 

to be not more than 10* of the maximum 
signal ) . 

We first consider a special form of eqn. (la) where phase 
differences are taken as zero i.e., < f > 3 = 0* We have 
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generated the results for 8j = 1, r 2, = 3 in figure Ala. The 

FFT of this curve has been carried out and given in figure Alb. la- 
this fugure the position of the spikes give o> t and the ordinate 
gives aj. For example if a^ is taken as unity a 2 and a 3 are 2 and 
3. This way it is that the signal is completely recovered. 

Various filter are used to eliminate z and b and highlight 
the importance of Uj f <*> 2 , Now we can select any frequency aiuL, 

using proper filters we can eliminate rest of them i.e., using 
band pass filters we can go for o> 2 , by low pass filter we can go 
for Wj and high pass filter for 

Now we will see effect of various windows in the regneration 
of signal : Here we have found that -rmet angular windows are not 
regenerating signals properly and with Hanning windows results are 
satisfactory. (Any window can't remove all the errors and noise so 
we will have to compromise in regeneration between amplitude and 
frequency. Since in our analysis frequencies have much important 
role to pay so we will select windows according to acceptable 
frequency regeneration). 

In the examples presented in Fig. la shows the signal having 

noise, dc component and three frequencies • which are 0.05, 1 and 7 

/ 

Hz. So we will use low pass filter of order 11 for elimination of 
7Hx frequency. Fig. Ala shows FFT magnitude response of crude 
signal which has three peaks corresponding to 0.05, 1 and 7 Hz 

frequencies. Fig. Ai c shows filtered FFT magnitude of signal with 
no peak for 7Hx frequency. 


APPENDIX 2 


SPECTRAL REPRESENTATION CF A TIME SIGNAL 

Discrete time signal x(n) is periodic if for some positive 
value of N, where N is period of the signal, 

x(n) = x (n+N) (2.1) 

The Complex exponential e *(2*/N)n ig periodic with period 

i • 6 * * 

*k <"> = *k+N (n) 

where 

= e ik ^ 2n,/ ^ J1 E (2.2) 

we write x(n) as x = £ e ik ^ 2,c / N ) n (2.3) 

n , . Ic 

k=o 


where the values of the function at different n values are 


N-l 

l *, 

k=o 


&T 1 2Kik/n 

Lt ® 

k=o * 


N-l 


N-l 


E « 

k=o r 


2*ik(N-l)/n 


If we wish to represent a Periodic Signal by a Fourier series, in 
the following we use eq. 2.3 
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I 


V 1 a e ik < 2,t / N > n = 1 + e ik(2*/N)n + ik(2n/N)n 

k=o 

(2.4) 

= 1 + a + ot 2 + . . .+ <x N-1 
A S 

where 

a = e ik(2Jr/N)n -.■* 

Above S can take on discrete values as 


S 


f H 





when a = 1 
when a * 1 


i 1c / / 1^[ \ 

Now a will be equal to 1 when e ' - 1 i.e. when k is a 

nultiple of N (k = 0, ± N, ±2N, In other words. 


I'* e ik(ZH/N)n = . 

n=o 

N 

0 

when k = 0, ± N, ±2N 

otherwise 

(2.6) 

Now * = I"' a k e ik < 2 ’' /N)/n 

k=o K 



- r 1 * ir(2*/w)/n _ 

" n 

n = o 

n 1 

C »-* 

5" 1 ik(2n/N)(k-r) 

t- k 

k=o J 

(2.7) 

N~1 

= L 

k=o 

N-l 

E 

n=o 

i (2 jc/N) (k-r )n 
e 
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N 


if (k-r) = 0, ± N,... 


N-l 



1° 


otherwise 


V _ fx J\ f since k = 0, 1,.... 

L, * iv "H i = o, 



where, 


a 


r 


l 

N 


N-l 

l 


n=o 


x 


n 


-ir(2x/N)n 


The Fourier series coefficients a^ are often referred to as 
the spectral coeff ici ent* These coefficients specify a 
decompose t ion of x r into a sua of N harmonically related complex 
exponentials, we if take k in the range from 0 to (N-l), we have : 


x n = Vo {n] + a l*l U) + • • • + a N _ 1 ♦ N _ 1 [n] (2.9) 

Similarly if we take k to range from 1 to N, we obtain 


X n = a i*i + a 2*2 + * * * " +rS N *N 


( 2 . 10 ) 


where 


Thus , 


* k In) a e 


ik(2«/N)n 


i(2K/N)n(k+N) = e i(2nKn/N) e i(2«nk) _ ^ [n] 


- * 

(J.U) 

On comparing equation (2.9) and (2.10), we get (using the fact 
= $ k+N InJ hence $ 0 [n] = 4> k In]) : 

*o = *N 
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that * k (n) 


i 


1 ftr 1 y i c 1 1 i k r&n^c over &ny set of N consecutive 

integers and using ^ = ♦j c+N *-We can conclude that 

a k = a k+N 


That is, if we consider more than N sequential values of k, the 
values of a k will repeat periodically with period N. In 
particular, since there are only N distinct complex exponentials 
that are periodic will period N, the discrete time Fourier series 
representation is a finite series with N terms r—Theref ore if we 
fix the N consecutive values of k over which we define the Fourier 
series given below, we will obtain a set of exactly N Fourier 
coefficients from equation (2.13) as follows : 


x 


n 


£ ‘ a e ik(2«/N)/n 


( 2 . 12 ) 


where , 

, - 1 „ -ik(2*/N)/n 

n ~ N n 

n~o 


(2.13) 


Above equation (2.12) is synthesis equation whereas (2.13) is 

analysis equation. Furthermore since the $ k [n] 

repeat periodically with period N as we vary k, so must the a^ 
vary (since a^ = a fc+N ) . 

Example 2.1 Consider the Signal 

x = Sin r n (2.1.1) 

n ° .in 

There are three different situations, depending on whether r 

o 

is an integer, a ratio of integers, or an irrational number , x ^ ^i s 


e ik(2n/N)n 
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in the third 


periodic in the .first two cases but not 
Consequently, the Fourier series representation of this signal 
applies only in the first two cases. 


2 Jl o v |, 

Case 1 is an integer N, * — = N 4 r Q = ~ 


4 x{nl = Sin r Q n = Sin ^ 


In 1 
2i 


i (2n/N)n -i(2n/N)n 

t; ** C 


On comparing it with 


N-l 


x „ = E a i 


ik(2n/N)/n 


n 


k=o 


4 a 


1 


2 i * 


-1 


1_ 

2i 


( 2 . 1 . 2 ) 


The remaining coefficients over the interval of summation are 
zero. These coefficients repeat in the period N. 


m 


Example 2.2 : Let x(n| = Sin (~-)n 


Period of repetition N=5. 


2* 


Case 2: — is a ratio of integers 
o 


✓ 4 IE 

_ N ^ 

*Z mmm mm r* 

J2«m| 

r 

m o 

l N J 

o 



* 

n 

* Star n 
o 

= Sin 


' inja 
l N 


3 - 


i_ 
2 i 


i ( 2«m/N )n -i<2ic*/N)n 
e — e 


N-l 


Compare with x = Y ’""'a' e 
n „ i 

m~o 


i*(2*/N)/n 


fiO 



a 


2 i ’ *- 


1__ 

2i 


Remaining coefficients over one period of length N are tero. 
The coefficients repeat with period N. 


2 JC 

Example : Let x|nj = Sin (-jr-)n 


■* a. 


'» a -3 


1_ 

2i 


‘3 " 2 i 
Period of repetition = 5. 


( 2 . 2 . 1 ) 


Exaup 1 e 2.3 

Consider, the discrete time periodic square were 



r 

0 

-N/2 

< n < 

" N 1 


From the figure x = * 

*1 

1 

“ N 1 

< n < 

N 1 

(2.3.1) 

-•*0 

V 


< n < 

N/2 


Use the relation a ~ 

f* 

I 

H 

N-l 

i x e 

n=o 

-ir(2x/N)/n 


* 





N, 


N 


l i. «- 


ir(2x/N)/n 


n=o 


i 

N 


e + e +ir(2*/H)/K£ e -^Fr(2n/N) + 


+ir(2*/H)/N, 


which is geometric series with (2Nj +l)_terms 


. ir<2x/N)N. 

a = ~ e * 

r N 


si 



1 e *( 2 / N ) N i e -ir(2K/N)(2N+l/2) p -ir(2Jt/N)(2N+l/2) 

«r , a r = N e -ir(2n/N)l/2 {e ir(2n/2N) _ e ir(2Ji/2N) 



v, ■ *r«**/N) * 0 

ff e 


(2.3.2) 


i.e. if r * 0 

if r = 0 then the sun is 



if r = 0 


Also a 


r+N 


This -expression for the Fourier series coefficients can be 
written more completely if we express the coefficients as samples 
of an envelope 


Sin ( (2Nj+i)r/2j 
N *k = Sinfr/2) 


r = 2«k/N 


(2.16) 


Discrete Fourier Transform Pair 


\ (■*) 

s 

N-l 

I 

n=o 

y d 

<»t r7*n**M* 

(2.17) 

Uj (nT) a 

1 

i 

H-l 

m=o 

-Hd 

f * 1 -i(2*n/N)/m 

{ NT J 


where. 







T is sampling integral 



Both w d and W d are periodic with a periond N. In particular 


W d. ( NT ) W d ( NT + T K ) 


(2.18) 


w. (nT) 

u 


= « . [ (n ± kN) T } 


Example 2.4 


2n _ N_ 2n 

r h1 o N 

o 


Also let 


2 * 


<^ d (nT) = Sin — na ’ { a ' not necessairly an integer) 


“d < nT > 


2 i 


i(2*a*/n)N -i(2Ha'n/N 

© **" © 


} 


-i (2Hna'/N) -i(2*n/N)m 

© © 


virs?.) 


r I j -i(2nn/N)(a'+a) -i(2nn/N )(«+■’ A 

1 U 1 * e J 


N-l 

l 

n=o 


i(2nn/N)(a'+a) -i(2nn/N )(»+■’ 


' i i •*-' 


i ) 2*/N ) (a * +a)N j i _ e " 1 < 2 */ N > m ’ + “ )N 


2i 


u _ c -i(2*/N)(a-»') u _ e -i(2K/N)(«a-m’) J 


Case 1 : Let a “ a’ and 2* be saall 


i(2*n/N)(a+a* ) s (1 _ p ) (j et ) where |p| is saall 


I - e 


•i <2*a/RHa+a 


' )j p is complex in general 
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TABLE - A2 


Some Examples of Filters 

Low Pass Filters - They permit lower frequencies to pass through 
and the high frequencies are stopped, with a transition zone of 
passband and stopband (of frequencies). 

We design such filters by choosing a finite, symmetric set of 
coefficients we select a digital filter of the form 

kN - — ■' 

y n = ? C k U n-K (C k = ^k* 

^ IN 

U n = U(n) sampled value. 

High Pass Filters - It passes the high frequencies and stops low. 
It can be called a difference between all pass filters and low 
pass filters. 

We design it by simply replacing f in L.P. filters by (1/2-f) 
it changes sign oLadd kt i so ni eg . 

Bandpass and Bandstop Filters - They permit/restrict a particular 
range of frequencies. They have two transition zones. 

/ 

NOTCH FILTERS - A very narrow band stop filter is called notch 
filter. Among other uses, a notch filter is used to remove the 
ever present 50 Hz that comes from electrical power distribution 
system. 

Consequently the first zero occurs where n=M. This function 
is such that it decreases rapidly with n and is negligible beyond 
its first zero n=M. If M is made very large, that is, the 
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integration segments are very small, then it will be the higher 
order coefficients (i.e .,3 large n) that are in the vicinity of 
the first zero. Let it be further assumed that the magnitude of 
these higher order coefficients is negligibly small. Then an 
approximation for the coefficients -- - 

M -i« tm 

D n " AT {1} £ f e n (2.2) 

m 

On the assumption that it is good for n < m and c n w 0 for n * M. 

Since there is no point in evaluating the coefficient for n > 
M-l, then the approximation for the Fourier series coefficient 
(now called Discrete Fourier Transform}. im—4*ken as 

. N-l +i« tm . N-l +i2*nm/N 

f m = f(t m> * T ^ D n e T ^ D n e (2 * 3) 

“ 1 n=0 n 1 n=0 n 

N-l +i« tm N-l -i2*nm/N 

D-~r-D(o ) * AT £ f e n = AT £ f e (2.4) 

n n n=0 " n=0 n 

where both m and n range from 0 to N-l. There are the definitions 
of the discrete Fourier Transform. Here the exponentials don't 
contain dimensional quantities, only the integers n,m,N. In this 
transform, both the time and frequency domains are discretized, 
and, as a consequence, the transform behav«s*periodical ly in both 
domains . 

PROPERTIES OF DISCRETE TRANSFORM 

t = t = mAT (2.5) 

HI 

0) = u n = 2nn/NAT (2.6) 
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also f = f(t ) is the input signal for the discrete Fourier 

■ B 

transform, while f(t) is the input signal for the continuous 
transform thus 

C ♦♦ D 

n 

That is, they both have the same magnitudes. This was 
achieved by choosing the appropriate scale factors of eqns. for 
approximate of the Fourier coefficient. 

FAST FOURIER TRANSFORM 

Final step in the numerical implementation is the development 
of an efficient algorithm for performing the -eummat ions of the 
discrete fourier transform. The fast fourier transform (FFT) is 
simply a very efficient scheme for computing the discrete Fourier 
transform. It is not the different transform. 

ALIASING 

The basic difference between the continuous and the discrete 
fourier transform is sampling (both rate and duration). An actual 
singla of very long duration must b« . tmmcated to a reasonable 
size before the discrete fourier transform is evaluated. 

The transform of the signal will therefore be dependent on 
the actual function F(t) and the way it is sampled. Note that if 
the sampling rate is every AT, then the highest detectable 
frequency using the FFT is 1/(2AT). Thus of AT is too large the 
highest detectable frequency using FFT is 1/2AT. Thus, if AT is 
too large the high frequency appears as a lower frequency cits 
alias) as shown in Fig. This phenomena of aliasing can be avoided 
if the sampling rate is high enough i.e., rate commensurate with 
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the highest significant frequency in the signal. 

From the practical point of view, this is also avoided by 
recording with the analog filters set so as to remove si gn*-$ leant 
frequencies above Nyquist. 

LEAKAGE 

The estimated spectrum is given only at discrete frequencies 
and this can cause another problem known as leakage. For example, 
suppose there is a special peak at 117 KHz and the sampling is 
such that the spectrum values are only at every 15 KHz. Then the 
energy associated with the 117 KHz peak leaks into the 
neighbouring frequencies, thus distoring their spectral estimates. 

Leakage is most significant in instances where the signal has- 
sharp spectral peaks and is least significant in cases where 
signal is a flat, bread - banded signal. In the analysis of the 
impact of structures, the signals generated generally don't 
exhibit spectral peaks (as in case of vibrations). However, it can 
become a problem if there are many reflections present in the 
signal . 

,P ADDING •*- 

Since discrete fourier analysis represets a finite sample of 
an infinite signal on a finite period. Then schemes for increasing 
the apparent period must be used. The simplest means of doing this 
is by padding. 

There are various ways of padding the signal and the one most 
commonly used is that of simply adding zeroes. 

The main consequence of padding is that since a large time 
window is used them higher resolution in the frequency domain in 


69 



achieved. It also makes the signal non-periodic over a large time 
domain. There effects can be seen. 

As regards padding with zeroes it should be pointed out that 
it really doesn's matter if the zeroes are before the recorded 
signal or after thus even if the signal never decrease to zero, 
padding is justified because it represents the only quiescent part 
of the signal. 
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Program to transfer time data from AD-3525 FFT Analyser^ 


REM 
REM 

REM ^ = = = = = = = == = =:=: = = = =:=: = =: = =: = =: = =: = = =:=:=: = = = =: = = = = = = = = = = = = = « = = = : 
REM 

CL S : LOCATE 4 , 7, 1, 0, 31 

REM PLAY w o2116e~9n0e“d-d-n0e-n0e-d-d-d-n0e-9n0e-d-d-n0e~” 
PRINT " PROGRAM TO TRANSFER TIME 
PRINT 
REM . 

CLEAR 
1 B IN 
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DATA 


. d eel . baa 
, 60000! 

IT1 « 60000! 


90 

100 

110 

120 

130 

140 

150 

160 

170 

ISO 

200 

210 

220 

225 

230 

240 

250 

260 

270 

280 

290 

340 

341 

343 

344 

345 
350 
355 
360 
375 
380 
385 
395 
415 
420 
425 
427 
430 

434 

435 
440 
445 
450 
455 
460 
462 
465 
470 


ittiw III = ouuuu { : IBINIT2 = IBINIT1 + 3: BLOAD ”bib.m’\ IBINIT1 

CALL IBINIT1 C IBFIND, I BTRG , IBCLR , IBPCT, IBSIC, IBLOC, IBPPC, IBBNA, 

I BONL , IBRSC, IBSRE, IBRSV, IBPAD, IBSAD , IBIST, IBDMA, 

I BEOS , IBTMO, IBEOT, IBRDF , IBURTF , IBTRAP ) 

CALL IBINIT2 C IBGTS , IBCAC, IBUAIT, IBPOKE, IBURT , IBURTA , IBCMD, IBCMDA, 

IBRD, IBRDA , IBSTOP, IBRPP, IBRSP, IBDIAG, IBXTRC, IBRDI, 

I BURT I , IBRDI A , IBURTIA, IBSTA*, IBERR% , IBCNT%) 

REM . ... assigning a unique identifier to device and storing it in DVM% 
BDNAMES = "KIK" 

CALL IBFIND(BDNAME$ , DVM* ) 

REM 

REM ....checking for error on IBFIND call. 

IF DVM% < 0 THEN GOSUB 820 
REM 

REM ....clearing device. 

CALL IBCLRCDVM*) 

IF IBSTA* < 0 THEN GOSUB 830 
REM 

REM ....setting sensitivity range to +1 OdBV . 

URT$ = "FRQ13 SNA8 SNB8 TDM1 TLN6 TGSO " : CALL IBURT(DVM*, URT$) 

IF IBSTA* < 0 THEN GOSUB 840 

URTS * "VDS1 VUL1 VAB1 VULO VABO TPVO.OOO TRPO TRG1 n : CALL IBURT( DVM-s ,URT$ ) 
IF IBSTA* < 0 THEN GOSUB 850 

PRINT "IBSTA* after initialising = IBSTA* 

CNT* = 4096 
CNT% * CNT% / 2 
DIM IARU5Q00): FLAG = 0 
PLAY "MBo0L32A~L6 4CL1 6BL6 4A + " 

INPUT "Press a key after collecting data H ;JUNK$ 

REM 

PLAY "MBoOL3 2A- L6 4CL1 6BL6 4 A> ” „ 

PRINT: INPUT "File name for Channel A data , r la$ 

PRINT: INPUT "File name for Channel B data --> ” , FLB$ 

CLS • FOR J = 0 TO CNT2*: IAR*(J) = 0: NEXT J 

URT$ = "THL1 TAR 00-1023? ” : CALL IBURT ( DVIU , URT$ ) 

IF IBSTA* < 0 THEN GOSUB 860 

F$=FLA4 

GOSUB 900 

REM 

PLAY ”MBoOL32A-L6 4 CL1 6BL64A+” 

CLS: FOR J = 0 TO CNT2S: IAR%(J) = 0: NEXT J 

URT$ = ”THL1 TBR 00-1023? ” : CALL IBURT ( DVM4 , URT$) 

IF I BSTA% < 0 THEN GOSUB 870 
F$=FLB$ 

GOSUB 900 
REM 

jjgfl Releasing FFT Analyser 

CALL IBCLRCDVM*) 

IF IBSTS* < 0 THEN GOSUB 875 
Vi = 1: CALL IBONLCDVM*. Vi) 

IF IBSTA ! : < 0 THEN GOSUB 880 
PRINT: PRINT 

PLAY ”MBoOL32A- L6 4CL16BL64A+” 

LOCATE 16, 21, 1, 0, 31: INPUT 
IF YESS = ”Y° THEN GOTO 10 

FNP 


DO YOU UISH TO CONTINUE CY/N) 


YESS 
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492 REM 
495 REM 

500 REM //////////////////Error Messages/////////////////// 
505 REM 

PRINT ” IBFIND call error”: RETURN 
PRINT ” IBCLR call error”: RETURN 
PRINT "1st I BURT call error”: RETURN 
PRINT "2nd IBURT call error”: RETURN 
PRINT "Channel A's IBURT call error”: RETURN 
Channel B’s IBURT call error": RETURN 


PRINT 
P 


PRINT ” IBRD I call error”: RETURN 


DATA. PLEASE UAIT 


820 
830 
840 
850 
860 

g 7 Q ******* *.**«**.*,%** _ w _ _ 

875 P RINT "IBCLR call error while releasing FFT”: RETURN 

87 6 ------ -- “ 

880 REN 
885 REN 

900 REN //////////////Subroutine for processing & writing data/////////////// 
910 REN 
915 REN 

920 CALL IBRD I (DVM4 , IARS(O), CNT4) 

923 IF IBSTA4 < 0 THEN GOSUB 876 
925 PRINT "NO. OF BYTES AFTER IBRD I = ”;IBCNT4 
930 PLAY "MBo0L32 A-L6 4CL1 6BL6 4A+” 

940 CLS : LOCATE 12,20,0,0,31: PRINT ” TRANSFERRING 
990 OPEN F$ FOR OUTPUT AS #1 
1010 FOR J = 0 TO (CNT24 -2) STEP 2 
1020 BYTE14 = IAR f .(J) AND 255 
1030 BYTE24 = (IAR4(J) AND 32512) \ 256 

1040 REM . . ^ 

1050 REM ...the leftmost byte, BYTE24, is the Sign-byte... 

1060 IF IAR4(J) < 0 THEN BYTE24 = BYTE24 + 128 
1080 BYTE34 = IAR4(J + 1) AND 255 
1090 BYTE44 = (IAR4(J + 1) AND 32512) \ 256 
1100 R EN 

1110 REM .the leftmost byte, BYTE44, is the sign-byte... 

1120 IF IARUJ + 1) < 0 THEN BYTE4 * = BYTE44 + 128 
1150 S% = (BYTE14 AND 128) \ 128 
1160 N14 = BYTE14 AND 127 
1170 L74 = (BYTE24 AND 128) \ 128 
1180 EE1 = Nil * 2 
1190 EE2 = L74 
1200 EE = EE1 + EE2 
1220 N24 = BYTE24 AND 127 
1230 FI = N24 / 128 
1240 F2 = BYTE34 / 32768! 

1250 F3 = BYTE 4 4 / 8388608! 

1260 F = FI + F2 + F3 
1270 IF (S4 = 0 AND EE = 0 
1280 VOLT = 0 
1290 GOTO 1350 
1300 EX = EE 

1310 TT = 2 * (EX - 127) 

1330 VOLT = (11 + F) * TT 
1340 IF S4 <> 0 THEN VOLT 
1350 REM 

1400 PRINT #1, VOLT 
1410 SOUND 200, .025 

1420 NEXT J 
1430 CLOSE #1 

1440 PLAY "MBT16001L8CDEDCDL4ECC" 

1450 CLS : LOCATE 11, 22, 1, 0, 31: 

1460 RETURN 


AND F = 0!) THEN 1280 ELSE 1300 


-VOLT 


PRINT ” DATA TRANSFERRED TO”, FS 



